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A procedure is described for the calculation of activity and osmotic coefficients which is hased upon a 
knowledge of the equilibria in solution and assumed single-ion activity coefficients. The procedure permits one 
to introduce chemical equilibria of various types (ion-pairing, complexation, hydration, and hydrolysis) into a 
model which can be used to calculate values of the excess Gibbs energy and the activity and osmotic 
coefficients. Both the Debye-Hiickel theory and Pitzer's expression are used to calculate the electrostatic 
contribution to the single-ion activity coefficients. Calculations have been performed on aqueous sulfuric acid, 
acetic acid, hydrofluoric acid, cadmium chloride, copper sulfate, and sodium carbonate. Properties which have 
been calculated are the excess Gibbs energy, the osmotic coefficient, the mean ionic activity coefficient, and 
Frank's single-ion activity coefficient function. Agreement between calculated and measured properties has 
been obtained up to molalities of ~ 1.0 mol kg~' 
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1. Introduction 

Equilibrium models have been used [1 — 15] 1 both for 
the prediction and for the correlation of activity and 
osmotic coefficient data in aqueous electrolyte solu- 
tions. These equilibrium models are particularly appro- 
priate when one is dealing with solutions which exhibit 
association, complexation, hydration, or hydrolysis. 
When applied to such solutions, they are superior to the 
use of a model that assumes the electrolyte in solution to 
be a fully dissociated strong electrolyte. A variety of 
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approaches has been used in these equilibrium models 
for treating these various types of equilibria in solution 
and for the calculation of the activity and osmotic coef- 
ficients. Several different types of functions for the cal- 
culation of the electrostatic contribution to the activity 
coefficient of the ions in solution have also been used. 
The purposes of this paper are to 1) describe a pro- 
cedure for (a) the calculation of activity and osmotic 
coefficients in aqueous solutions that uses a generalized 
approach for treating the equilibria in solution, and (b) 
the calculation of the Gibbs energy properties, 2) clarify 
the distinction between the stoichiometric and species 
quantities, frequently a source of confusion in the litera- 
ture, and 3) explore the effects of parameter variations in 
the model on calculated values of the thermodynamic 
properties as applied to several representative types of 
electrolyte solutions. Using the equilibrium model, one 
can also calculate the values of the activity coefficients 
of individual ions. This permits one to then calculate 
values for Frank's [16] single-ion activity coefficient 
function. 
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2. The Model 

A fundamental idea used here is the distinction be- 
tween the stoichiometric components of which a system 
is composed and the particle constituents, or species, 
which are introduced to account for the properties of a 
solution. Stoichiometric quantities will be designated by 
"st" when not otherwise clearly indicated. Species or 
particle quantities wiJJ always be designated by a " " " 
placed over the quantity. For example, for the solution 
formed from 0.1 moles of pure H 2 SCXr(0 and 1 kg/M] 
moles of H 2 0(^), the stoichiometric molality of H 2 S0 4 
is m<H 2 SO 4 )=0.1 mol kg -1 , m(SO! _ )=0.1 mol kg -1 , 



and m(H*) = 0.2 mol kg -1 if one views the electrolyte 
as being completely dissociated. However, if one con- 
siders the equilibrium: SO+~"(aq)-rH~(aq)=HSO! _ (aq), 
adopts a K of 99" for it [17], and uses a Debye-Hiickel 
type expression (see eq (5)) for the activity coefficients 
of the ions with a B parameter equal to zero, one calcu- 
lates m(H + ) = 0.134 mol kg -1 , m(SO3")=0,G342 mol 
kg- 1 , m(HSO 4) =0.0658 mol kg" 1 , and m(H 3 SO 4 )=0.0 



! Throughout this paper the activity will always have units of mol 
kg -1 . The equilibrium constants, formed as products and quotients of 
activities, will have units of mol kg -1 , or kg inol~',or, for symmetrical 
reactions, will be dhnensionless. For sake of brevity, the unils for 
equilibrium constants will be omitted. 
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mol kg -1 . Physical quantities other than amounts of sub- 
stance can also be viewed as particle quantities. In par- 
ticular it is the distinction between G" and & c * which 
forms the basis of the model presented herein. 

The general system to be considered is formed from 1 
kg/Mi moles of water and n^(/=2 to JV ) moles of 
other components. In terms of a particle model, the 
system can be viewed as being formed from h\ moles of 
water and n L (i=2 to iV s ) moles of particles or species. 
The amount of water (»]) will always be designated by a 
subscript "1." The components and the species other 
than water will be designated by subscripts "f" and "i," 
respectively, with i and />2. 

Each component £ in the solution is represented as 
C z ' c A " where C and A are the reference cations and 
anions, respectively. The charges of the cations and an- 
ions are z^ c and z^ a , respectively; the ion-numbers are 
v fc and v fa< respectively. It is important to note that the 
choice of reference species is arbitrary; e.g., for a solu- 
tion of aqueous sulfuric acid, the reference species could 
be selected as either two H + ions and one SO 2 , - ion or as 
one H + ion and one HSO4 ion. In the former case, 
sulfuric acid would be represented as (H+^CSO 2 .^! and 
in the latter case as (H^ifHSO^)!, While the choice of 
reference species is arbitrary, it will not affect the 
amounts of the various species one calculates given the 
same set of equilibrium constants and allowed spe- 
cies. However, as will be seen later, the choice does 
affect the values of several of the stoichiometric thermo- 
dynamic properties. 

The Gibbs energy of the solution in the stoichiometric 
representation is given by 



G = (mfkg)G,+ ^Jv ec G fc +v ea (3^)mt 



(1) 



where m f is the molality of water in pure water and G (c 
and Gfz are, respectively, the chemical potentials of the 
reference cation and anion of the f' h component in the 
solution. In terms of an equilibrium or species model of 
the solution, the Gibbs energy is given by 



:2 



G=.2«iGi 



(2) 



If the equilibrium model is accurate, the Gibbs energy of 
the solution calculated using eqs (1) and (2) will be 
identical and the chemical potentials and the activities of 
the \' h species will be the same in both the stoichiometric 
and species representations. The development to follow 
will start with a description of the solution in terms of 
the equilibria assumed to be present in solution, an as- 
sumed expression(s) for the activity coefficients of the 
solute species in solution, and the calculation of the 



activity of the water using the Gibbs-Duhem equation. 
It is important to note that there may be several equi- 
librium models of a solution which yield agreement with 
calculated properties. Thus agreement between calcu- 
lated and measured properties does not, in the absence 
of direct molecular information, prove the correctness 
of the model used. 

The equilibria in solution are described by a series of 
chemical equations: 



JVfc 

2*jk Sjk = 0, 



j=i 



k = 1 tO iV e 



(3) 



where iSj k is the j'* species in the k rt equilibrium, N k is the 
number of species in the k"" equilibrium, JV e is the total 
number of equilibria, and f jk is the number of moles of 
species S jk participating in a given equilibrium; ? jk is pos- 
itive if 5j k is a product and is negative if 5j k is a reactant. 
The equilibrium constants are: 



JVk 



K k =TT a? 



", jk 



j=i 



k= 1 to N c 



(4) 



where a >k is the activity of species S >k . Since Sj is equal to 
the product of the molality of the i' A species {m^ and its 
activity coefficient (fi), the complete formulation of the 
equations which describe the equilibria in solution re- 
quires that some assumptions(s) be made concerning the 
form of the %'s (i>2) in solution. In this paper two 
different expressions for the ys's will be used: 



fn% = -A m m m /(\+BI m ) 



(5) 



and 



€n %= -%A$ tn /(\+bT n ) + Q./b)t ■nil+bV' 2 )] (6) 

Equation (5) is the Debye-Hiickel equation with an ex- 
cluded volume or "ion-size" parameter B in the denom- 
inator. Equation (6) is the leading term of Pitzer's ex- 
pression for fi [18]; he has set b equal to 1.2. A m and A$ 
are Debye-Hiickel constants, where A m — 3A$= 1.17642 
kg 1/2 mol" 1/2 at 298.15 K. Equations (5) and (6) can be 
extended by the addition of the expression: 



JV S N t N s N % N, 

2 2 Xifhj+ 2 2 2 HijkmjWk 

]=j j=l i=] j=t k=l 



(7) 



where \jj and juty k are, respectively, the interaction pa- 
rameters for pairs and triplets of particles. We shall later 
return to the subject of the extensions of eqs (5) and (6) 
and to several other aspects of the choice of an expres- 
sion for fj. It should be noted that / is calculated as 
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/=(l/2) 2 rhizf 

i = 2 



(8) 



We have formulated the equilibrium equations using 
an extent of reaction variable £. Thus, the amount of the 
i' A species in solution is given by 



«i = (m f kg)6 M + 2 (v ^" ^S.v + v ^ a n ^8,>) 






i>l 



(9) 



where £ k is the extent of reaction variable for the k" 1 
equilibrium and S,j is the Kronecker delta; r is an integer 
which serves to identify the reference cation and anion 
of each of the components in the solution. The first 
summation term on the right side of eq (9) specifies the 
amount of species i which is formed in the absence of 
any equilibria in solution; the second summation term 
specifies the contributions, which may be positive or 
negative, to the amount of species i in solution from the 
equilibria in solution. If the water is a participant in the 
equilibria in solution, «] must also be calculated using the 
above equation. Hydration numbers can be introduced 
directly using 



«i=«,(eq (9))- 2 h-fr 

1=2 



(10) 



where h t (eq(9)) is the amount of water calculated using 
eq (9) and ft-, is the number of waters of hydration at- 
tached to the i'* species. If the water participates in the 
equilibria in solution or if hydration numbers are intro- 
duced, «] kg -1 will not be equal to m f. The ih\, are given 
by 



m i =« i (mf/n 1 ), i>2 

The activities of the species are given by 
ai=m$ lt 1>2 



(11) 



(12) 



Thus in a model where hydration is introduced, both m, 
and a { will be affected by changes in both n, and «; (i>2). 
To obtain a numerical solution of eq (4) , it is neces- 
sary to make some initial guess for the activity of water 
if it is a participant in the equilibria in solution; we have 
generally used a value of unity. Thus, having formulated 
the simultaneous nonlinear eq (4), which neces- 
sarily include eqs (8) to (12) and eq (5) or (6), one is left 
with a numerical problem to obtain a self-consistent 
solution of these equations. It is assumed that, while 
such a solution may be difficult to obtain for large sys- 
tems, a unique solution does exist and that one now has 



values for /, m, and «j for i>l and for y-, and a x for i>2. 
The activity of the water can now be calculated by 
application of the Gibbs-Duhem equation stated in 
terms of excess properties of the species 



i=2 



Use of eqs (U) and (13) leads to 



i=2 



(13) 



(14) 



It is also necessary to adopt some conventions concern- 
ing the limits of a t and % The conventions used herein 

are ay— ►l and, for i>2, %— »1 as 2 /Mr- »-0; also, y? (i>2) 

i=2 

is defined to be equal to unity 3 . These conventions, to- 
gether with the definitions of the activity a„ the activity 
coefficient y„ and the definition of the excess Gibbs 
energy of the i' h species given in the following equations 



a 1 =exp((G i -Gf)//ir] 



y i =a i /m i 



,i>l 



,i>2 



(15) 



(16) 



and 



Gr=Gi-G; d ,i>i (i7) 

lead to 

G?=RTShy, (18) 

The introduction of eqs (5) and (18) into eq (14) yields 
&S=(2A m RT/m tB 3 )[(l +BI m ) 

-2<?«(H-5/ 1/2 )-(l+5? 1/z >-'] (19) 

Similarly, one obtains, using eq (6) instead of eq (5), the 
relationship 



GT=(2^ R T/m f ) r? 3/ 7(l + bT /2 )] 



(20) 



We now consider an ideal reacting solution in which 
yi (i>2) is equal to unity at all temperatures, pressures, 

3 This is the convention most frequently used in the description of 
aqueous electrolyte solutions. It is based upon the molality scale. A 
different convention based upon the mole fraction scale is commonly 
used for the description of non-electrolyte solutions. 



254 



and compositions and which also allows for the pres- 
ence of both equilibria and hydration in solution. Appli- 
cation of the Gibbs-Duhem equation to such a solution 
leads to 



of component / is given by 



ye±={y^"-y f : fi ) 



•v.yAv 



(26) 



in h?= ■ 



- 2 rh-Jml 

1=2 



(21) 



and this in turn with eq (15) leads to 



G? = G° x -(RT/mf) 2 Oil 

i=2 



(22) 



Since an excess property is defined as the difference 
between the real and the ideal, it follows from eqs (15) 
and (22) that 



'na { = &f/R r - ( 2 iWi/m f ) 

i = 2 



(23) 



If the equilibrium model is an accurate representation of 
the solution, a, is equal to a { . We thus have a procedure 
for the computation of the activity of the water which 
starts with an equilibrium model of the solution and an 
assumed expression for the activity coefficients in the 
solution. Since it is necessary to make an initial guess as 
to the value of the activity of the water, the calculation 
should be repeated using the value of oj from the pre- 
vious iteration until convergence to within a given toler- 
ance in its value is obtained. The (stoichiometric) os- 
motic coefficient is calculated as 



<J>= 



-(mf/2 v^mAina. 



(24) 



As was the case for the osmotic coefficient, the value of 
the mean ionic activity coefficient is dependent upon the 
reference species selected for a given component. 

Other stoichiometric Gibbs energy properties can be 
calculated in addition to the activity of the water, the 
osmotic coefficient, and the mean ionic activity coeffi- 
cient of the /'* component. Additional properties of the 
water are calculated as follows: 



Gf=RTina t + (RT/mf) i v e m e 

f-2 



G x = G1+RTina\ 



and 



y!=ai/^i 



<fc, =(/««,)/( /BAD 



(27) 



(28) 



(29) 



(30) 



Note that eqs (29) and (30) are definitions of the activity 
coefficient of the water (yi), and of the rational osmotic 
coefficient (4> x ), respectively. These two quantities have 
not been frequently used in the literature. The properties 
of the solutes are: 



where m/ is the stoichiometric molality of component 
i and Vf is equal to (v^ c + i-vJ. Note that v f is unity for 
a non-electrolytic component and that, as stated earlier, 
the value of the osmotic coefficient is dependent upon 
the choice of the reference species selected for each 
component in the solution. 

The stoichiometric activity coefficients can be calcu- 
lated using the principle that the chemical potential is 
independent of any representation of it. Equating a„ to 
a„ the stoichiometric single-ion activity coefficient is 
given by 



yt=m$/m\ 



(25) 



where mf is the total stoichiometric molality of the i'* 
species in solution. The mean ionic activity coefficient 



and 



G?=v ( RTiny e± 



y,=exp(G$*/RT) 



a e =m f y ( 



G f -G}=RTina e 



a ( ±=m {± y i± 



(31) 
(32) 
(33) 
(34) 

(35) 



where the mean ionic molality of component i is de- 
fined as 



mt±=i(y e jne) Vt '^feJney n Y Nt 



(36) 
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Finally, the total properties of the solution are 



and 



G"= 2 n,G? 



G-G°= 2 n/Qf-G?) 



(37) 



(38) 



Equation (27) is the stoichiometric analogue of eq (23). 
Equations (29) and (30) are, respectively, the definitions 
of the activity coefficient of the water and of the rational 
osmotic coefficient. The steps used in the overall com- 
putational procedure are summarized in figure 1. 

It is interesting to consider the consequences of the 
ideal behavior of all of the solute particles in the species 
representation, i.e., %= 1 for i>2. Application of eqs 
(21), (23), and (24) to such a solution leads to 



AT, ^ N c 

cf>= 2 fhj 2 v .m? 

i=2 f=2 ' 



(39) 



the case when there are no chemical interactions (i.e., 
association, hydration, or hydrolysis) present in solution 
and the stoichiometric reference species chosen are the 
only ones present in solution. 

An additional interesting feature of this model is that 
properties of individual ions and species are calculated 
(eq (25)). Unlike the stoichiometric properties given in 
eq (24) and eqs (26) to (35), there are presently no experi- 
mental data available with which one can compare these 
calculated values. Frank [16] has defined a quantity S ± 
characteristic of single-ion properties for a binary elec- 
trolyte: 



8 ± =(yyy;y 



(41) 



Application of eqs (25) and (26) leads to 



The definition is easily extended to multicomponent sys- 
tems as was done for the case of the mean ionic activity 
coefficient (see eq (26)) and values of S^± can also be 
calculated if the equilibria in solution are known. 

An alternative way of viewing this model uses the 
definition of the excess Gibbs energy 



y ■ e ± = {m c /m3'«(fkjmy 



(40) 



G tx =G-G li 



(42) 



If the stoichiometric and species representations are 

identical, then <f> and y f± will be equal to unity. This is and its analogue for the species representation 



Figure 1 — Steps for the calculation of the amounts of species in solu- 
tion and of the stoichiometric Gibbs energy properties of the solu- 
tion. 



10. 



12. Calculate the other stoichiometric 
properties, eqs (27) to (35). 

11. Calculate «f> using eq (24) and y ± 
using eqs (25) and (26). 

Make a choice of stoichiometric reference 
species. 



9. If Oi is equal to the initial guess of Si, go to step 10; 
if not, return to step 6. 

8. Calculate a { =a h eq (23). 
7. Calculate Gf, using eq (19) or (20). 

6. Solve eqs (4) in a manner self-consistent with eqs (8) to (12) and with 
eq (5) or (6). 

5. Make an initial guess for a\. 
4. Choose an expression for y a eq (5) or (6). 
3. Specify the values of K K , eq (4). 
2. Specify the stoichiometric amounts of each component in the solution. 
1. Set up the equilibrium eqs (1). 
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G !ex =G--G id 

Since the Gibbs energy is independent of representation 



(43) of these or any other choices for y u it does impose one 
important constraint on such a choice, namely that 



G~=&* + &«-&< (44) or equivalents, 



Similar equations hold for G" 
tion of ideality, 



<f > 1 . From the defini- 



&?=&\+RTfnmx 



(45) 



Equations analogous to the foregoing expression and to 
eq (39) exist for the stoichiometric components of a 
solution. Thus, introduction of eqs (22) and (45) and 
their stoichiometric analogues into eq (44) yields 

G e *={2 «iGf} 



i=l 



•>V S 



+{« 1 [G?- J Rr(.2m/mf)] 






+ j£Jv /e nv(G°c+*r/iiv, e iM,) 



(46) 



Note that G? = G° in the above equation. 

The three terms in { } on the right side of eq (46) 
correspond, respectively, to G e \ & d , and G ld . Inspection 
of the terms for G" 1 shows that the differences between 
these quantities involve two factors: 1) a difference in 
Gibbs energies, i.e., G? terms multiplied by n- K and «i, and 
2) entropic terms, i.e., the {RTItiifnm^ terms for the 
solute particles and the (it 7*2/71;) terms for the solvent. 
If one views G™ as the electrostatic or ionic contribution 
to G", the stoichiometric excess Gibbs energy is seen to 
also consist of energetic and entropic contributions 
which are formally accounted for with this model. 
While eq (46) could be used to compute G™ directly, it 
is numerically preferable to use the computational 
scheme outlined earlier (see fig. 1). 

In this paper, two different expressions have been 
used for y, (see eqs (5) and (6)). While classical thermo- 
dynamics has little to say about the correctness of either 



(d Sn y/afl| J -)r^«^ i =(3 /» 7A«drj,,n m 



(47) 



(48) 



Thus, while it is tempting to try to assign a different 
value olBorb to each species in a solution in eqs (5) and 
(6), respectively, to do so would violate eqs (47) and 
(48). However, the extension of eqs (5) and (6) using eq 
(7) does not violate this thermodynamic constraint. 
Note that the use of eqs (5) or (6) for 7, does not allow 
for the introduction of specific-ion effects attributable to 
long-range electrostatic interactions. These effects can 
be introduced by the use of eq (7). Specific-ion effects 
attributable to chemical equilibria are accounted forin 
the equilibrium part of the model. Expressions other 
than those in eqs (5) and (6) could be used to represent 
the electrostatic part of y„ 

The long range electrostatic contributions to the 
Gibbs energy properties are introduced via eq (5) or (6) 
and the use of an equilibrium modeJ. Other interactions 
accounted for in this model include: 1) chemical inter- 
actions, 2) hydration, and 3) volume exclusion effects. 
The attractive chemical interactions are accounted for 
by the use of the equilibrium constants for the processes 
which describe the equilibria in solution. These pro- 
cesses can involve ion-pairing, complexation, and hy- 
drolysis. The effects of hydration are accounted for 
either by the introduction of equilibrium constants for 
specific reactions involving hydration or by the use of 
hydration numbers for each species in solution. The use 
of hydration numbers reduces the value of h\ in eq (10) 
which in turn has consequent effects on the m^, 7, % a m 
and other properties. Volume exclusion effects are rep- 
resented by the B or b parameter in eqs (5) or (6). Short 
range repulsive forces between particles can also be ac- 
counted for using the Xy and/or u^ parameters. It is 
worth noting that there are similarities in the effects that 
changes in certain parameters have on thermodynamic 
properties. Specifically, an increased value of B (or b) is 
similar to the introduction of a positive Xy or ^ and also 
to the introduction of hydration effects. Physically this 
should be the case since the excluded volume for a hy- 
drated ion is larger than for one that is not hydrated, 
Also, a negative X^ is similar to an association between 
particles i and j. The remainder of this paper will discuss 
the application of this model to several aqueous salt 
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solutions containing representative types of chemical 
interactions. 

3. Results and Discussion 

"We now compare the results of calculations using this 
model with experimental data and examine the results of 
perturbing the various input parameters in the model, 
i.e., the single-ion activity coefficient expressions, the 
assumed equilibrium constants, and the assumed state of 
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Figure 2— ClockwisE from above, calculated activity coefficients, os- 
motic coefficients, and excess Gibbs energies of aqueous HiSO* at 
298. 15 K as a function of molality. Ka is varied: 99 for the solid line, 
89 for the dashed line, an 109 for the dotted line. Equation (6), 
Pitzer's expression for % was used in calculating each of the 
curves. The ions were considered to be unhydrated. 



hydration of the species in solution. To do this, calcu- 
lations have been made on the following aqueous elec- 
trolyte solutions: sulfuric acid, acetic acid, hydrofluoric 
acid, cadmium chloride, copper sulfate, and sodium car- 
bonate. 

The results of calculations on aqueous sulfuric acid 
are shown in figures 2 through 4, where the chemical 
equilibrium considered is 



H + (aq) + SOi-(aq) = HSO^Caq) 
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Figure 3 — Comparison of calculated osmotic coefficients for aqueous 
H 2 S0 4 at 298. 1 5 K. The squares are the correlated values of Pitzer, 
Roy, and Silvester [10]. The solid line was calculated using a value 
of K A = 99 and Pitzer's expression for % The dashed and dotted 
lines were calculated using the same value of K^ and values of B 
equal to zero and 2.5 respectively, in a Debye-Hiickel expression 
for 'Ji (eq (5)). The ions were considered to be unhydrated. 



The value of K A is 99 at 298. 15 K [17]. It is seen that the 
calculated values of y±, tf>, and G ex are relatively insen- 
sitive to moderate variations in K A up to molalities of 0.2 
mol kg -1 . However, as seen in figure 3, significant 
changes in the osmotic coefficient are produced by per- 
turbing the value of the B parameter. Use of a value of 
B equal to 2.5 in a Debye-Hiickel expression for y { (see 
eq (5)) produces good agreement with the experimental 
osmotic coefficients (the squares in fig. 2d) up to a mo- 
lality of ^0.2 mol kg -1 . The effects of variations in 
hydration numbers is shown in figure 4. The minimum in 
the osmotic coefficient curve cannot be produced by 
variation in the B parameter, but, as seen in figure 4, a 
minimum is observed when hydration is introduced. 

The equilibrium considered in the description of 
aqueous acetic acid (HAc) is 




molality/ (mol kg 



Figure 4— Calculated osmotic coefficients of aqueous H2SO4 at 298.15 
K in which the extent of hydration of the ions is varied. All three 
curves were calculated using a value of 99 for Kj and a Debye- 
Hiickel expression for ^j with B equal to 2.0. The dotted line was 
calculated assuming that all of the ions were not hydrated. The 
solid and dashed line were calculated, respectively, by assuming all 
three ions (H + , SO4 ~, and HSO4) to be hydrated with three and 
four waters each. The squares are from the correlation of Pitzer, 
Roy, and Silvester [10]. 



H + (aq)+F-(aq)=HF a (aq) JC = 1.44XI0 J (C) 



and 

HF (aq) + F-(aq) = HF7(aq) ,*T=2.71 (D) 

For aqueous cadmium chloride the stepwise equilibria 
considered are: 

Cd 2+ (aq)+Cl-(aq)-CdCl + (aq) ,tf = 85 (E) 

CdCl + (aq) + Cl-(aq)=CdCl?(aq) ,#=2.71 (F) 

CdClS(aq)+Cl-(aq)=CdCll(aq) ,^=0.53 (G) 



H + (aq)+Ac-(aq)=HAc°(aq) ,* = 5.96xl0 4 (B) and 



For aqueous hydrofluoric acid, the equilibria are 



CdCl7(aq)+Cl-(aq)=CdCir(aq) ,^=4.3X10-* (H) 
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The equilibria considered for aqueous copper sulfate and 



are: 



Cu I+ (aq)+SOj-(aq)=CuSO?(aq) ^=250 (I) 

CuSOS(aq) + Cu 2+ (aq)=Cu I SO^(aq) ^=5 (J) 

CuSOS(aq) + SOl-(aq)=Cu(SO^-(aq) ,K=5 (K) 

For aqueous sodium carbonate the equilibria considered 
are: 

COi-(aq)+H 2 0(/) 
=HCOl(aq)+OH-(aq) ,K= 2.16X1Q - * (L) 
HCOl(aq)=CO?(aq)+OH-(aq) ^ = 2.34 X 10" a (M) 



3-d 
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Figure 5— The osmotic coefficients of aqueous acetic acid at 273, 15 K_ 
The squares are the experimental osmotic coefficients reported by 
Harris, Thompson, and Wood [21]. The solid curve was calculated 
using a value of K B = 5.96x 10* which was obtained from the val- 
ues of AG', LH°, and AC^ for process (B) given by Larson and 
Hepler [19]. The solid curve was also calculated using either Pit- 
zer's ft or a Debye-Huckel ft where B varied from zero to 2.0. In 
this figure and in all subsequent ones, the ions were considered to 
be unhydrated. 



H 2 0(0=H + (aq)+OH-(aq) JT=1.0X10- 



(N) 



AH of the above equilibrium constants refer to a tem- 
perature of 298.15 K. with the exception of the value for, 
process (B) which refers to 273.15 K. The Gibbs ener- 
gies of formation given in the NBS Tables of Chemical 
Thermodynamic Properties [17] were used to calculate 
the above values for equilibria (C), (D), (K), (L), and 
(M). The value for process (B) was calculated from the 
data for acetic acid tabulated by Larson and Hepler [19]; 
the values for processes (E), (F), (G), and (H) are those 
given by Reilly and Stokes [12]; and the values for pro- 
cesses (I), (J), and (K) are those given by Pitzer [9]. The 
Debye-Hiickel constants recommended by Clarke and 
Glew [20] were used in all calculations, 

The results of the calculations are shown in figures 5 
through 9. For acetic acid, since the species ionic 
strength (1) is very low (it has a value of 0.044 mol kg" l 
at P'= 1.0 mol kg -1 ), the choice of the expression for y-, 
makes very little difference. Near agreement with the 
experimental osmotic coefficients is obtained to a 
molality of ;=0.2 mol kg -1 . The difference between the 
measured osmotic coefficients and the calculated ones 




molality/ (mol kg 



Figure 6— The osmotic coefficients of aqueous HF at 298.15 K. Values 
of K c = 1442 and Kq = 2. 63 were used to calculate the solid curve 
together with Pitzer's expression for ft. The squares are from the 
correlation of Harrier and Wu [22]. 
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can be attributed to the formation of dimers and trimers 
of acetic acid [21] and cannot be explained either by the 
introduction of hydration or by the use of different B 
parameters in eq (5) or by the choice of a different 
expression for j\. 

For aqueous HF, CdCl 2 , CuS0 4 , and Na 2 C0 3 , agree- 
ment of calculated with measured properties can be ob- 
tained to molalities between 0.6 and 1.0 mol kg ' by the 
variation of the B parameter in eq (5) and, for the case of 
CdCl 2 , using only eq (6) for y,. Neither the use of eq (5) 
nor eq (6) is able to produce the minima in the calculated 
values of <{> or y ± which is observed for many electro- 
lyte solutions. These minima can be produced, however, 
either by the introduction of hydration or by the use of 
eq (7) to extend the equations for % It should be noted 
that for Na 2 C0 3 solutions the osmotic coefficient does 
not approach the usual limit of unity as m st approaches 
zero mol kg '; instead it approaches a value of 1.395 
[14]. This is a consequence of eq (39) and the presence of 
equilibria (L), (M), and (N). 

Calculated values of 5 ± are shown in figure 10. The 
fact that values of S ± for acetic acid and for copper 
sulfate are essentially unity is attributable to a very near 
cancellation of terms in eqs (25) and (36). It should be 
noted that the value of 8 ± for Na 2 C0 3 , unlike the other 




0.4 0.6 
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Figure 7-The activity coefficients of aqueous CdCl 2 at 298.15 K. The 
successive formation constants given by Reilly and Stokes [12] 
were used in doing the calculations. The solid line was calculated 
using Pitzer's expression for yi and the dashed line was calculated 
using the Debye-Hiickel y, with B set equal to 2.0. The squares are 
based on the measurements of Reilly and Stoke [12], 



systems shown in figure 10, does not approach a value of 
unity as the molality approaches zero mol kg -1 ; the 
minimum value of S ± for Na 2 C0 3 occurs at a molality of 
0.0040 mol kg -1 . While there are presently no experi- 
mental values of 8 ± available with which to compare 
our calculated values, this property is potentially mea- 
sureable [24]. 

In summary, an equilibrium model for aqueous solu- 
tions has several important applications: 1) the Gibbs 
energy properties can be reliably estimated at low 
molalities if the appropriate equilibrium constants are 
known, 2) an equation of state can be generated which 
is appropriate for a particular type of solution, 3) 
amounts of species in a given solution can be calculated, 
4) single-ion activities can be calculated, and 5) as was 
done here, effects of variations in the equilibria, state of 
hydration, and electrostatic contributions to the Gibbs 
energy properties can be investigated. A natural exten- 
sion of this model is the calculation of enthalplies, heat 
capacities, and volumes of aqueous solutions. 



The author thanks Drs. Graham Morrison, Ralph L. 
Nuttall, and, in particular, Robert H. Wood, for useful 
and stimulating discussions on the topics of this paper. 
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Figure 8 — The osmotic coefficients of aqueous CUSO4 at 298.15 K. 
The squares are the experimental data of Miller et al [23]. The three 
curves were obtained using values of Ke = 250 and Kp=Kq = 5. 
The solid line was calculated using Pitzer's expression for % the 
dashed and dotted lines were obtained using the Debye-Hiickel yi 
with B set equal to 2.0 and 5.0, respectively. 
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Figure 9— The osmotic coefficients of aqueous Na 2 COj at 289.15 K. 
The squares are from the correlation of Vanderzee [14]. The solid 
line was calculated using Pitzer's expression for ?i and values of 
&Lt Km< and Kfj from reference [17]. The dashed line was calcu- 
lated using these same values of A!" and a Debye-Huckel ■?, with B 
equal to 1.0. 
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